Modeling DNA conformational transformations on the mesoscopic scales 
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The approach for the description of the DNA conformational transformations on the mesoscopic 
scales in the frame of the double helix is presented. Due to consideration of the joint motions of 
DNA structural elements along the conformational pathways the models for different transforma- 
tions may be constructed in the unifying two-component form. One component of the model is 
the degree of freedom of the elastic rod and another component - the effective coordinate of the 
conformational transformation. The internal and external model components are interrelated, as 
it is characteristic for the DNA structure organization. It is shown that the kinetic energy of the 
conformational transformation of heterogeneous DNA may be put in homogeneous form. In the 
frame of the developed approach the static excitations of the DNA structure under the transitions 
O ■ between the stable states are found for internal and external components. The comparison of the 

' data obtained with the experiment on intrinsic DNA deformability shows good qualitative agree- 

ment. The conclusion is made that the found excitations in the DNA structure may be classificated 
^— ^ . as the static conformational solitons. 

PACS numbers: 87.15.-v, 87.14.Gg, 46.25.-y, 63.90.+t 

I. INTRODUCTION 

The study of the conformational transformations of DNA double helix is the line of understanding the mechanisms 
of the living systems functioning. Last decade particular attention attracts the conformational transformations of 
I ' DNA macromolecule, which include the displacements of the structural elements in the frame of double helix with 
O" 1 relatively large deviations from equilibrium positions. The amplitudes of these displacements are not small, but on the 
other hand - not so large to destroy the double helix organization. Such DNA transformations in the frame of double 
\ helix we will name as mesoscopic, leaving the term "large-scale" for the displacements of the structural elements that 
J> ■ can destroy the helix. The DNA mesoscopic transformations show themselves in the formation of intermediate A-B 
conformations in the complex of DNA TATA-box with proteins [1-4], the motifs of A-form in ligand-bound DNA 
(see [5] and cited therein), the junctions of B- and A- forms in DNA fragments bound to the protein or drug [6-9], 
the bistable flexibility of A-tract [10-12], the base-pair preopening in oligonucleotide duplex [13], the overstretching 
of DNA helix under high tension [14,15], and in a number of unusual DNA conformations [16]. In mentioned above 
' conformation transformations the displacements of DNA structural elements are accompanied with the deformation 
f*^; , of DNA macromolecule as a whole, but without the destruction of the double helix. 

The investigation of the mesoscopic transformations is of great interest because these DNA rearrangements directly 
connect with the regulation of the genes activity [1,17]. From a complex structure of DNA macromolecule the 
■ experimental study of its conformational mobility meets with some obstacles [18]. One of the productive way to 
understand the conformational resources of the double helix is the use the phenomenological approach to modeling 
its conformational transformations. The exploring such models as "helix-coil" [19-21] or "elastic rod" [22-26] speeds 
up and extends the investigations of the conformational properties of the nucleic acids, and allow to interpret a set of 
experimental data on DNA melting, helix bending and torsion, loop formation and higher order folding. The model 
approach has definite advantages comparing to various all atoms computational methods. From the large number 
of atoms in a macromolecule the use of the computational approach is restricted by the fragments of DNA helix 
not longer than 20 base pairs [27,28]. At the same time under successful model construction the phenomenological 
approach allows to describe the structural transformations in macromolecules of the real size and with the account of 
surrounding factors. 

For modeling DNA mesoscopic transformations it would be advantageous to use an experience in the modeling 
of large-amplitude mobility of the structural elements in DNA double helix, and the study of the possibility of the 
localized excitations realization [29] . On this direction quite a number of models have been constructed (see review 
[30], collected articles [31] and book [32]). But not all of these studies achieved real results because of the significant 
simplifications of the models and the difficulties of tacking into account such principle properties of DNA structure, 
as nucleotide content and sequence. 

Recently, a set of the models for the reduced description of the internal mobility of DNA structural elements on 
the large scales has appeared [33-37]. The authors of this studies divide the atoms of the double helix on the atomic 
groups - the new "structural" elements of the system, and try to find the potential field which could describe the 
interactions between them. This is a very complicated task which cannot be resolved for DNA as a whole. 
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In the present work the another way is proposed. In modeling DNA conformational transformations on the meso- 
scopic scales the four-mass model approach is used. This approach was developed earlier [38-40] for the description 
of DNA double helix conformational vibrations. The achieved agreement between theory and experiment in the 
description of the DNA low-frequency spectra [39,40] show the accordance of the four-mass model to the described 
mobility of the double helix. It would appear reasonable that the displacements of DNA structural elements with 
small amplitudes are the precedes of the mesoscopic and large amplitude transformations. It is taken into account 
also that under double helix transformations the structural elements move jointly along the conformational pathways. 
Therefore, starting from the four-mass model for DNA internal mobility the models for definite pathways of the double 
helix transformations are constructed. Such an approach allows us to reduce significantly the degrees of freedom of a 
system and to choose the adequate potential functions for described force field. 

Using the proposed approach the models for the description of DNA conformation transition of A-B type and double 
helix opening arc constructed. It is shown that different conformational rearrangements of the structural elements in 
the double helix may be described by the unifying two-component model. The accounting of the interrelation between 
the model components is discussed as an important part of the description of DNA transformations. In the frame 
of the two-component model the static excitations are found for two types of the conformational transformations of 
the macromolecule: transition from the ground state to metastable and transition between the equivalent states in 
the conditions of DNA bistasbility (for example, A-B equilibrium). The obtained results are compared with with the 
experimental data on DNA deformability. 

II. MODEL CONSTRUCTION 

The DNA molecule consists of two strands which are formed by the linked phosphate groups and deoxyribose (or 
sugar) rings. To the sugar rings the nucleic bases are attached. The DNA bases have different content of atoms, that is 
why the primary structure of a macromolecule is heterogeneous and may be used for writing the genetic information. 
In the natural conditions DNA strands form the double helix - DNA secondary structure, where the nucleic bases of 
different strands are connected by hydrogen bonds in the complementary pairs A-T and G-C. 

The DNA rearrangements on non-small scales in the frame of the double helix have a character of transitions from 
one stable state to another [18]. Under conformational transformations the disposition of the structural elements 
in double helix monomer units changes, and the macromolecule as a whole suffers definite deformation. The form 
of macromolecule deformation and the type of the conformational state of the double helix are observed in the 
experimental studies [1-16]. So, for modeling the mesoscopic transformations of DNA double helix it is important that 
the constructed model will include the components for the description of the internal conformational transformation 
and the external deformation of macromolecule. For construction the models of definite DNA transformations let us 
consider the four-mass model approach. 

A. Generalized four-mass model 

According to approach of four-mass model for the DNA monomer link [38-40] the basic conformational degrees of 
freedom, which characterize the mobility of the structural elements of a double helix, are the following. First of all, 
the displacements of the nucleosides (nucleic bases together with sugar rings) as rigid pendulums " suspended" to the 
backbone (Fig. la). Secondly, intranucleoside transformations: the displacements of nucleic bases under the changes 
in the sugar rings conformation. The third type is the displacements of a nucleotide as a whole (nucleic base + sugar 
+ phosphate group). 

For the description of DNA structural dynamics in the approach [38-40] two masses of the pendulum-nucleosides 
(m) and two masses of the backbone atom groups (m ) (Fig. 1) were used. The masses of nucleosides have the different 
values according to the kind of the nucleic base. The backbone masses mo include the group of atoms: PO2 + 20 
+ C5/ and are identical along the chain. The necessary structure parameters of the model, such as the length of the 
physical pendulum-nucleosides - Zo, and the angles of their position in the pair plane - 60 (Fig. 1) calculate from the 
known double helix conformation. For DNA B-form the calculated parameters of the model are shown in the Tabic I. 
As seen, four-mass model takes into account an essential property of the DNA macromolecule - heterogeneity of its 
monomer content. Really, the model masses and geometrical parameters of the pendulum-nucleosides depend on the 
nucleic bases contained in its constitution. 

Let us consider the possibility of using the four-mass model for description of the mesoscopic transformations of 
DNA double helix. We will write the expression for the energy of generalized four-mass model without envisaging the 



smallness of displacements: 

E= ^{ m oRL + m 4 r2r i + [/(r 4 ,„;R 4 ,„)}. (1) 

In the expression (1) the summation is made over all monomer links (n) of the double chain (i — 1,2). The radius- 
vectors Rj i7 j(i) and r i n (t) describe the displacements of the mass centers of the backbone groups and nucleosides, 
respectively, in each chain. The R and r denote the time derivatives. The function U is the potential energy of 
conformational transformations of a system. 

Taking into account the structural organization of a macromolecule, the potential energy of the n-th monomer link 
would be presented as: 

U(n) = U {1) {n) + C/ (2) (n,n - 1) + C/ (3) (n; n ± 1). (2) 

The first term in the sum (2) describes the energy related to the internal rearrangements in n-th monomer of the 
macromolecule as such. In the four-mass model this energy may be written as: 

= [/ 1 (1) [R 1 -r 1 ]+C/ 2 (1) [R 2 -r 2 ] + 
+ ^ 1) [r 1 -r 2 ] + (7| 1) [R 1 -R 2 ], 

where and are the energies of the nucleoside displacement with the respect to the backbone group in each 

strand, is the energy of the H-bond stretching in the base pairs, and is the energy of the interaction 

of the backbone groups of different helix chains. But, at once the term may be omitted because interaction 
between backbone groups in a monomer link is much less than the same interaction along macromolecule chain (the 
corresponding distances differ more than twice [18]). The interactions along the backbone chains will be taken into 
account in another potential terms. 

The term in the sum (2) describes the energy of the interaction along the macromolecule chain. This energy 
term may be written in the form corresponding to the traditional approach of the elastic rod model with the accounting 
interactions between the nearest neighbors. The view of the term would be concretized for definite coordinates 
of the conformational rearrangements. 

The term U^ 3 ' describes the energy of the interrelation between the intermonomer transformations (conformational 
changes), and the changes of the configuration of the helix as a whole (deformation of the macromolecular chain). 
The term is important under the consideration of the mesoscopic displacements of the macromolecule structural 
elements, because DNA structure is relatively soft, and under the change of the monomer configuration the position 
of monomer link in a macromolecule chain changes also [18]. 

The energy expression (1,2) is sufficiently complicated for the direct use. Firstly, the vector character of the 
displacements which requires the construction of the corresponding force field for the macromolecule. Secondly, the 
heterogeneity of the energy terms, that makes the analytical study practically impossible. 

But, as known [18,25,27,41-47] under rearrangements of DNA structure in the frame of the double helix the displace- 
ments of the structural elements of the monomer link are interdependent, and DNA structural elements move jointly 
along the definite conformational pathways. Let us take this fact into account and consider the model constructing 
for the basic conformational transformations of DNA double helix: the conformational transition of A-B type and 
the helix opening. 



B. Conformational transitions 



Under modeling DNA conformational transitions let us introduce the displacements of the center of masses of the 
pair of nucleosides (r p ) and the center of masses of the backbone groups (r q ), as well as the displacements of the 
nucleosides with respect to each other in the monomer link (5), and the same for the backbone masses (A ): 

niiri + m 2 r 2 . 
r P = m p — > m P = mi+m 2 ; 

(3) 

8 = ri - r 2 , n P = mi • m 2 / m p , 



where m p is the mass of pair of the nucleosides in the monomer link, and u, p is the reduced mass of the nucleosides. 
Accordingly, for the backbone masses we have: 

r q = (Ri + R 2 ) , m q = 2m ; 

A = Ri - R 2 , n q = m /2 

In the designations (3) and (4) the expression for the n-th item of the kinetic energy of a system takes the form: 

1 ( • 2 • 2 ^ 

K n = - |m p r p 2 + m q v q + u, p S + /j, q A j . (5) 

For the further advance it is useful to put out the displacements of the mass center of a monomer link (R) and 
displacements of the mass center of a pair with respect to the backbone masses (r): 

R = "V> + "V\ M = m 1 +m 2 + 2m ; 

r = r p — v q , to = nip ■ m q / M . 

Here M and to are the mass and reduced mass of a monomer link, respectively. Taking into account (6) an expression 
for the kinetic energy takes the following form: 

K = \Y,{ M ^n +mi n+^l+V q K)- (7) 

n 

Notice, the expression (7) can be simplified some more. At the description of the dynamics of conformational 
transitions not all terms of the kinetic energy are equally important. As was mentioned above, the potential energy 
associated with the last term in the kinetic energy (7) may be omitted {U^ ~ 0). Thus, from the point of classical 
mechanics A is the cycle coordinate, and it can be ignored here. So, it is possible to neglect the last term in the 
expression (7). 

Let us take into account also, that in such structural transformations as conformational transitions between the 
forms of B- and A- families of a double helix, the nucleic bases in the pair displace together without significant mutual 
shifts. In any case, under conformational transition both in the initial and in the final helix forms the bases remain in a 
configuration of a complementary pair [5,18,41-43,47]. So, we will assume that under DNA conformational transitions 
it is true: 8 « 8 « 0. Therefore, an expression for the model kinetic energy for the description of conformational 
transitions can be written as the following: 

n 

The kinetic energy in the form (8) remains still sufficiently complicated for the modeling. And, as the last step 
in model construction, let us take into account that in the conformational transformations of the macromolcculc 
the masses M and to have the primary directions of the displacements. As it is known, under the conformational 
transitions of A-B type the mass center of DNA monomer link has the most displacement in the plane which is 
orthogonal to the helix axis [18,41-43]. Thus, for the degree of freedom of the monomer mass center it is convenient 
to pass to the displacements in the plane XOY (we put the helix axis along 0Z, Fig. 1). In this case we have: 

MR 2 = Mt 2 + /^ 2 , 

where T and ip are the displacement and the angle of rotation of the link in the plane XOY, / = I z is the moment 
of inertia of the monomer link with respect to a helix axis. 

The necessity of taking into account the displacement on T and torsion on ip simultaneously is determined by 
the concrete conditions of the transition. Both types of deformation are observed in the experiment studies [5,18]. 
Under definite conformational transitions these degrees of freedom change in a correlative character to minimize the 
distortion of base pair stacking. The including in the model both degrees of freedom do not cause the difficulties. For 
the purposes of the present study under description of A-B transition we will use one coordinate - T to denote this 
type of the transformation. 

In its turn it is also known [18,43,47,48], that under conformational transitions the base pairs displace transversely 
to the helix axis, along the axis OX (Fig. lb). Thus, for the component r it is expedient to pass to the displacement 
along the trajectory of the mass center motion: r — > u. 



(4) 



(6) 



In such a way an expression for the kinetic energy of the A-B transition in terms of T and u can be written in the 
scalar form: 

K tr = 1 -J2{MTl + mul}. (9) 

n 

Under constructing the kinetic energy of the system one fact has attracted the attention at once. After passing 
in the kinetic energy to the description of the joint motions (beginning from the expression (7)) the appeared mass 
coefficients don't depend on the heterogeneity of the monomer content of DNA macromolecule. The calculated values 
of the masses are shown in the Table II for A-T and GC pairs. 

The potential energy for the conformational transition in the selected coordinates consists of the energy of helix 
bending (external component) and the energy of base pair displacements (internal component). 

The energy term in expression (2) in this case comes to the energy of the pair displacements relative to the 
backbone. These displacements occur due to the changing in the sugar ring conformation [18,41,43,47,48], and the 
potential energy along the coordinate u corresponds to the transition from one sugar form to another. We will describe 
this energy by the function 

The term in expression (2) which describes the energy of the .ff-bonds stretching in the base pairs has to be 
omitted because the conformational transformations of A-B type, as was mentioned above occur without significant 
change of the distance between nucleosides mass centers (S « 0). 

The potential energy of the interaction along the macromolecule chain (term in expression (2)) includes the 
energy of the bending on the component T which is bound to have the view akin to the bending energy of the elastic 
rod, and the energy of the change of pairs stacking interactions (displacement u). 

Therefore, the expression for the potential energy for the description of the dynamics of the conformational transi- 
tions of DNA macromolecule takes the form: 

Utr = 2Ef T [ T " _T "- 1 ] 2 +* : »[' 1 » _ "»-lf 



+ $[u n } + W[u n ;r n±1 \y (10) 

In the expression (10) ky and fc„ are the elastic constants of the interactions along the macromolecular chain for 
bending and pair displacement components, respectively. 

The term with the function W describes the interrelation between the base pair displacement in n-th monomer and 
the corresponding displacement of the n-th monomer with respect to (n — l)-th and (n + l)-th monomers (the term 
U^ 3) in the expression (2)). Here: T„±i = [T„ + i - T„] + [T„ - T„_i]=[T n+ i - T„_i]. 

The energy expression (9,10) may be used for modeling the dynamics of DNA conformational transitions. It is 
important to note that owing to passing to the description of joint motions, the mass coefficients in kinetic energy 
became non-sensitive to the nucleotide heterogeneity in DNA. Moreover, the potential energy of the conformational 
transition dynamics became also non-sensitive to DNA nucleotide content. The effect of DNA heterogeneity may 
reveal only in some sensitivity of model elastic constants (k) to the nucleotide sequence in macromolecule because of 
known sensitivity to the sequence of the potential energy of interactions along the DNA chain [18]. These properties 
of the model are in common agreement with the known data about non-sensitivity of DNA A-B transitions to the 
nucleotide content [41,48]. 



C. Opening 

Other structural transformation that has attracted considerable interest in DNA physics is the pair opening in DNA 
double helix. This process offers the conformational transformation, as a result of which the hydrogen bonds in the 
complementary pairs disrupt, and the nucleic bases turn out from the double helix [18]. 

It is considered [27,46] that the most probable trajectory for double helix opening is the turning of bases of the pair 
around backbone chains to the sides of the helix grooves (so called " opening" , the direction of rotation on the angles 
0\ and 02, Fig. lb). The other possible pathway connected with the stretching of base pair along the H-bonds (OY 
axis on Fig. lb, "stretching") is assumed as less probable because of known strong rigidity of the backbone chains in 
the direction of OY axis [42] . The " stretching" process can realize under the temperature increase in the conditions of 
melting. The modeling of the base pair stretching in double helix was made in a number of works [30,49-52]) with the 
help of one-component model. In the present study the "opening" trajectory will be considered, as the most probable 
in the natural conditions where the mesoscopic transformations occur. 



Let us consider the form of kinetic energy for description of the double helix opening. In the four-mass model 
(1,2) we will introduce the displacements of mass centers of the nucleosides with respect to the backbone, and the 
displacements of mass centers of the nucleotides as a whole. For the displacements of masses in the n-th link of the 
z-th DNA strand we will write the following equations: 

= m, : r ? +m n R,, ^ = ^ + ^ . 

P l = Ri - r i5 ft = mi- mo/ Mj , 

where Mj and ft are total and reduced masses of the nucleotide in the z-th strand. 

Since the main process in the opening transformation is the extension the H-bonds in the base pair, it may be 
assumed that under helix opening the nucleosides remain rigid. We will consider that in DNA opening the mobility 
of the nucleosides is brought to the turnings as a whole around the backbone. Thus: 

ft 2 — k 2 e\\ (12) 

where k is the reduced length of the pendulum-nucleoside and 9i is the angle of its deviation from the equilibrium 
state (Fig. lb). 

Under the conditions (12) the kinetic energy of a system takes the form: 

K=^[ M ^ln+^ln]- (13) 

The expression (13) remains still complicated for the modeling of the DNA opening. Let us pass to the consideration 
of the joint motions of the structural elements in the double helix link. With this purpose we will single out the 
displacements of the mass center of both nucleotides (mass center of monomer link) and the relative displacements of 
nucleotides in the link. These displacements and the corresponding expressions for the masses have the form: 



R = 



= M iPi+ M 2E> 2 ) M = M 1+ M 2 ; 



(14) 



€ = P1-P2, M = Mi-M 2 /M. 

Let us introduce also the joint turnings of the nucleosides: 

n = JiOi-J202 } j = Jl + j 2 (Ji = M 2 ); 

(15) 

a = 6>i+0 2 , j = Ji-Ji/J. 

Here the parameters J and j have a meaning of the joint and relative inertia moments of the nucleosides in the pair. 
At writing (15) it was taken into account, that the angles Q\ and 62 are disdirective (Fig. lb). The angle a correspond 
to the known parameter " opening" in the classification of DNA transformations [53] . 

Taking into account the expressions (14) and (15) the kinetic energy of the opening can be presented as: 

K °p = I E { Mk « + + + 3°l } • (16) 

n 

The values of the parameters for the kinetic energy (16) are calculated and adduced in the Table III. It is significant 
that these parameters are practically independent from the kind of DNA pair. Thus, as in the case of the description 
of conformational transitions, both for helix opening the kinetic energy of the heterogeneous macromolecule assumes 
the homogeneous form after passing to the joint motions of the structural elements. 

Then, for further simplification we will take into account that in the opening process the motions of the structural 
elements take place along the definite pathway. We will presume the uniform distribution of energy among the moved 
nucleosides, and in the wake of [27,46] that the opening occur in such way that: Ox rts 2 and |f2| -C |<r|. So, we will 
suppose that: ft ~ 17 ~ 0. 

It is also known that when the base pair opens to the side of the helix groove, the monomer link as a whole displaces, 
and the macromolecule bends [44-46]. On this way the nucleotides in the pair displace jointly, that is corresponded 
to: p x k, p 2 , |£| <c |R|, and R — > Y, where Y is the coordinate of the displacement to the side of the helix grooves 
(Fig. lb). ^ 



Thus, an expression for kinetic energy of the double helix opening transformations acquire the two-component form: 



(17) 



The expression for potential energy for the opening pathway may be presented in the form which looks like the 
expression (10): 



U, 



op 



\ \- Yn ~ Yn - 1 \ 2 + k ° [ CT ™ - <7n-l] 



+ *[j n ]+W[(T n ;Y n±1 \j. (18) 

In the expression (18) ky and k a are the elastic constants of the interaction along the chain, $[cr] - the potential 
energy of the pair opening which also includes the nucleoside turning around the backbone chains, and the term W 
describes the energy of the interrelation between the components a and Y. 

The shape of the potential $[er] must allow for the description of transition from the ground state of the double 
helix to the metastable state. Considering the form of the potential <j>[cr] it must be taken into account that on the 
pathway of opening the interaction of the opened bases with surrounded water have lead to the formation a set of 
metastable states of the base pair with molecules of water. Such states (preopened) have been recently found in the 
quantum- mechanical calculations [54,55], and observed in the experiment [13]. These results show that the double 
helix opening is the stage- by- stage process, and the potential <j>[er] has bistable (or more correctly, multistable) form. 
In the case of the description of the transition from a close state of the double helix to the preopened one the potential 
<&[<r] must have the shape of double- well function. 

When modeling the potential energy of DNA opening it is important to allow the dependence of the energy of the 
pair opening on the kind of a pair (the difference in hydrogen bonding of A-T and G-C). In other words, in the case 
of opening the energy <&[<r] is sensitive to the base pair content in DNA. 



III. UNIFYING MODEL FOR DNA INTERNAL MOBILITY ON THE MESOSCOPIC SCALES 

As it is seen from the results of the previous section of this work, the relatively simple two-component models 
can be constructed for description of various DNA transformations in the frame of double-helical state. The models 
for different DNA mesoscopic rearrangements are very similar. They include the external component which is the 
degree of freedom of the elastic rod, and the internal component which is the characteristic coordinate of structural 
transformation. This is the basic construction that may be modificate by including the additional terms - components 
of both types. Let us write the expression for the energy of DNA mesoscopic transformation in the unifying form: 

E = ^Y^{ Mti n + k 1 [R n -R n - 1 ] 2 

+ mfl + k 2 [r ] 2 + <t>[r n }+W[r n ;R n±1 ]y (19) 

In the expression (19) R and r are the degrees of freedom of external and internal components, respectively, M and m 
- the full and reduced monomer masses, k\ and ki - the elastic constants of the interactions along the macromolecular 
chain for external and internal components. 

The function $[r] describes the potential energy of the internal subsystem on the pathway of a conformational 
transformation from one stable state to another. Hence, it must have the shape of the double well (Fig. 2). The form 
of the double well is usually known from the experimental data. For the conditions when one of the helix form is 
more stable (for example at the physiological conditions the most stable is the B-form of the double helix [18]) the 
function <J>[r] must be taken in the non-symmetric view with two non-equivalent stable states: ground and metastable 
(Fig. 2a). That is true both for the conformational transitions such as A—B, and for the transition in the preopened 
state. If the conformational transition takes place under equilibrium of the states (for example in the interval of A-B 
equilibrium), the function $[r] must have the symmetric view with two equivalent states (Fig. 2b). 

The term W reflects the energy of the interrelation between the components, and may be written in the form: 



W = X F[r n ]lR n+1 -R n - 1 



(20) 



Here \ is the coefBcicnt, the potential function F(r) describes the change in the interrelation energy under the 
transition of the internal component from one stable state to another. By its physical meaning the potential F(r) is 
the function that increases from the ground state of the system and decreases near the another stable state. In the 
case of the description of the transition between the equilibrium states this function must have the symmetric form 
in accordance to physical equivalence of the states. A sample view of this potential function is shown in Fig. 2a and 
2b.. 

The two-component models of type (20) were used for modeling the transformations of bistable chains of different 
nature. Firstly such model was used under study of the dynamical properties of the systems with H-bonds [56,57], 
where the bistable chains with equivalent stable states were considered. The two-component model of bistable macro- 
molecule with non-equivalent stable states was introduced for the study of the nonlinear dynamics of local A-B 
transitions in DNA [58-61]. The same models were used for the description of DNA melting [62], and the nonlinear 
dynamics of the endothcrmic transitions [63] . 

Let us pass to the continuum approximation, that is frequently used in analytic studies of models. In this case it is 
supposed that conformational excitation embrace some piece of the macromolecule chain which is much larger than 
the chain step. In the continuum approximation the expression for the energy (19) may be written in the form: 

M [r 2 + s\R' 2 ) + m (f 2 + sir' 2 ) 



J 2h 



+ $(r) + 2 X hF{r)R' 



(21) 



where the expression for the interrelation energy (20) was also taken into account. In the expression (21) R = R(z, t) 
and r = r(z, t), R' and r' are the derivations on z, s\ = kih 2 /M and s\ = fe/i 2 /m. 

Notice, when modeling DNA conformational transformations the sensitivity of the elastic constants to the nucleotide 
sequence was not considered. Really, the elastic constant ki (and to some extent fei) reflect the value of stacking 
interactions between the neighbor pairs in the helix. In its turn the stacking depends on the kind of the interacting 
pairs (see Rcfs. [18,27,48] and the references therein). This sensitivity of the elastic constants has an effect on the 
form of the fine structure of conformational excitation, and can be the object of special consideration under study of 
such effects on the quantitative level. 

The equations of motion for two-component excitation with the energy (21) have the form: 

R = s 2 i?" + Xi^V; (22) 
, „ 1 d<& dF . ,. 

where \i — X^/M and X2 = xh/m ~ the constants of the interrelations. 

The equations (22,23) describe the dynamics of the two-component lattice. It is very much the same that the dy- 
namics of one-dimensional molecular crystal [64] . The significant distinction is the including the terms with coefficient 
X- As was mentioned above, the interrelation between the components for such " soft" lattice as DNA macromolecule 
have been taken into account because of direct evidences of the experiment. 

For better understanding the role of interrelation let us consider the dynamics of two-component model (21) in the 
neighborhood of it ground state. For the small deviation from the ground state: r = r$ + £, £ <C ro and R ~ p <C 
the expressions for the potential functions may be written as: $ « ao£ 2 and F w &o£, where ao — \^pi\r a , and 
b = jfplro are the first nonzero terms in the potential expansions on £. 

The equations of motion (22,23) for the small deviations from the ground state have the form: 

P = s\p" + xib^'; (24) 

I = s%e - -a Q i - X 2b p' ■ (25) 
m 

We will find the decisions of the equations (24,25) in the form of oscillating waves: £ = £o expi(wt — kz) and 
p = poexpi(wt — kz). After the wave substitution the equations (24,25) lead to the dispersion equation which may 
be written as: 

(uj 2 - s 2 k 2 )(cu 2 - s 2 k 2 - 2a /m) = k 2 b 2 X iX2 ■ (26) 

For very small wave vector k (k 2 ~ 0) from (26) the two limit frequencies follow: uj 2 w lu 2 — 2a^/m and ui 2 w 0. 
Thus, in the two-component system (21) we have two type of phonon vibrations. One, with the limit frequency 



ujq = -^^pr\ ro - optical branch, that determines the dynamics of the internal component. Second type, with the 
u> — > under k — > 0, is the acoustic branch. 

The optical branch correspond to the definite optical vibration of the double helix, considered in the linear four- 
mass model [38-40] . That may be used for more precise selection of the actual degrees of freedom for modeling the 
mesoscopic mobility, and the evaluation of the parameters of the potential functions. 

For small k (k 3 and k 4 ~ 0, the quadratic approximation) the expression (26) is equivalent to the following: 

(cj 2 -c 2 k 2 )(u; 2 -c 2 k 2 -u; 2 )=0 7 (27) 

where c\ = s\- si and c| = s\ + si, si = 6^x1X2/^0- 

So, in the quadratic approximation for the dispersion law the analyzed system is the same that is true for usual 
two-component lattice [64] but with the force constants c\ and c 2 . The constants include the contribution of the 
energy of interrelation (s» ) . As seen from their values the interaction of the sublattices leads to the softening of the 
external sublattice - the degrees of freedom of the elastic rod, and to the greater rigidity of the internal sublattice. 
Notice that the effect of the component interrelation is larger than the frequency of the corresponding conformational 
vibration is lower. As known (see [38-40] and reference therein), the frequency of DNA conformational vibrations 
are sufficiently low (<100cm^ 1 ) So, the effect of the component interrelation is important, and must be taken into 
account under the modeling DNA mesoscopic transformations. 



IV. STATIC STATES OF THE TWO-COMPONENT MODEL 



In this section we will consider the static states of the unifying model (21), which may be compared with experiment 
to demonstrate the correctness of the approach. To the static states we will assign the stable static states of the system 
and the static excitations which can be stable under some conditions. The stable static excitations of the system 
may be observed as static deformations of the real macromolecule. So, let us determine the static states of the 
two-component model (21). 

For more general character we will study the static states of the system without definition of concrete expressions 
for the potential functions. In accordance with the shape of the potentials we will assumed that the potential energy 
of the system will have three extremums: two minimums (r , r 2 ) and one maximum (ri) between them. In the ground 
state r = ro (B-form for DNA) we will consider $(ro) = F(ro) = 0. In the metastable state (r = r 2 ) $(r 2 ) 7^ 0, 
F{ r 2) 7^ (Fig. 2a). For the conditions when the macromolecule have two equivalent stable states (Fig. 2b) the 
function F(r) for r 2 give F(r 2 ) = 0. 

For providing of the bistabile form of the model energy it is necessary to satisfy the definite conditions for the 
cxtrcmum points [65] : 

$(ri, 2 ) > j^F 2 (r h2 ) , (28) 

that is an inequality for r\ state always, and for r 2 state is inequality in the case of the non-equivalent states and - 
equality for the case of equivalent states. 



A. Static excitation in the macromolecule with metastable states 



Let us find the static excitations of the macromolecule with metastable states in the monomer link. We will 
assume: r = r(z) and R = R(z), and R'(z) = r/h. Here the value t(z) is coined for the determination of the 
macromolecule deformation on the external component (torsion or bending, in dependence of concrete model). The 
equations described the static excitations follow from the Eqs. (22,23) of the system: 

R" + Xify = 0; (29) 
1 rf$ „ dF „. „ 

In the Eqs. (28,29) xi = x/kih, X2 = x/foh and n = k 2 h 2 . 
After one time integration of the Eq. (29) we obtain: 



R' + XiF(r) = 0, 



(31) 



where the constant of integration put to zero according to the initial conditions for the stable state. The deformation 
of the external component of the macromolcculc chain may be written from the Eq. (31) as: 



(32) 



As seen from expression (32) the deformation of the chain in the model is proportional to the function F(r). 

After substitution the expression for R' in the Eq. (30) and one time integration we obtain the equation for the 
internal component: 



which may be rewritten in the form: 



r' 2 --<S>{r)+ Xl x 2 F 2 (r) =C, 



Q(r) = 0. 



(33) 



The equation (33) have the view of the energy of the mechanical system with the potential energy: Q(r) 
-^{r)+XiX2F 2 (r)-C. 
The solution of the Eq. (33) may be found as the conversion of the integral 



r(z) 



du 



r(0) 



(34) 



The view of the effective potential Q(r) is determined by the concrete form of the potential functions $ and F. 
When fulfilling the condition (28) the effective potential Q(r) have the form of a double hump (see Fig. 3a). 

Let us consider the possible solutions of the equation (34). We will find the excitations with a restricted trajectory 
and the asymptote of one of the stable states of the system because such excitations are the most stable in the real 
conditions. For the potential Q(r) with two non-equivalent humps we will search the excitation with the asymptote 
of the metastable state r 2 (Fig. 3a). In this case it is convenient to present the polynom Q(r) in the form: 



Q(r) 



l (r 2 - r) 2 (r - r?i)(r - r) 2 ). 



(35) 



In the expression (35) Q 2 ^ is the combination of the model parameters, r 2 , r\\ and r\ 2 (r]i < r\ 2 < r 2 ) are the zeros of 
the polynom Q(r) with the meaning C = — ^(r 2 ) + XiX 2 F 2 (r 2 ) for the boundary conditions of the metastable state 
(r = r 2 and r' 2 = 0). We will interest in the excitations with the trajectory restricted by the interval: r\ 2 < r < r 2 
(see Fig. 3a). 

For the polynom Q(r) in the form (35) the integral equation for the function r(z) have the form: 



r(z) 



du 



(r 2 - u)y/{u-r]{)(u-r]2) 



(36) 



The integral (36) have the table view. After integration we obtain the following relation: 

1 



Or, after rewriting: 



In ■ 



(r 2 -r)(r?2 - m) 



V( r 2 -Vi)(r-m) 



V ( r 2 - mjir - m) t = ±QmZ\Z(r 2 - m)( r 2 - m) ■ 



exp(<7 m z) (r 2 - r) (n 2 - ni) 



V (r 2 - vi)( r - m) - V ( r 2 - m)( r - vi) 



(37) 



where q m = ±Q m y/{r 2 - 772)^2 - m)- 



For convenience let us designate: r a — r 2 — rj 2 and r b = 772 — rji. The values r a and r b reflect the form of the 
potential Q(r) (Fig. 3a). Their substitution to the Eq. (37) lead to the expression: 

(r 2 - r) [r b exp(g m ) + 2r a + r b ] - 2r a (r a + r b ) 



= 2y/r a (r a + n) [(r 2 - r) 2 - (r 2 - r){2r a + r b ) + r a {r a + r b )} . (38) 

After raising to the second power of the whole expression (38) it is easy to obtain the form of the internal component 
of the static excitation: 

r(z)=r 2 - r cx (z) , (39) 

, x_ 2r a (r a + r b ) 

r b ch(q m z) + 2r a + r b ' 

As seen, at z — > (the middle point of the excitation) r cx — ► r a and r(z) — > -q 2 . At z — > 00, r ox — ► and r(z) — > r 2 
(Fig. 3b). 

The external component of the excitation according to the Eq. (31) have the form: 

r^(z) = -^F[r 2 -r ex (z)]. (40) 

At z — > 00, r cx — > 0, r cx — > r 2 and at z — ► 0, r cx — > r a , r ox — *■ r a = -^-F[r 2 - r a ], where > r a > t 2 . 

Thus, as seen from the Eqs. (39,40) the static excitation has the form of the bell with the width that is proportional 
to g" 1 ~ (r 2 . + r a r b )~2 and the height - r a . At the middle point of the excitation the components have the values 
that are close to the same for the ground state (Fig. 3b). On the edges of the excitation the components have the 
values of the metastable state. 

The obtained results may be used for comparison with the observed static deformation of the macromolcculc. The 
set of data on the static conformational deformation in DNA was published recently. So, in the work [66] it was found 
the structure of DNA-protein complex. The reconstruction of DNA structure in this complex [5,67] show that protein 
has induced B — > A transition in DNA, and the DNA fragment is bent and looked as the bell. Besides that, the 
conformation of the double helix in the center of deformation is close to ground state (i3-form), and on the edges - to 
A-form. 

Notice, that the static excitation (39,40) is commonly agree with the observed form of DNA deformation. As 
in experiment [5,66,67] the obtained static excitation has the maximum deformation on the edges of the fragment 
(r = T2), and the minimum - at the central part (r = r a , which t 2 < r a < 0) (Fig. 3b). It is also important, that the 
theory gives for the internal component the conformation of the metastable state on the edges of the excited fragment, 
and the conformation closed to the ground state for the central part. That is also in accordance with the experiment 
[5,66,67]. 



B. Static excitation in the macromolecule with equivalent states in the monomer link 

Let us consider the static excitations of the two-component chain (21) with equivalent states in the monomer link. 
As in previous case the form of the excitation for internal and external components is determined by the Eqs. (32-34). 
For the equivalent states in the monomer the effective potential Q(r) in the Eq. (33) has the shape of the double 
hump with equal humps (Fig. 4a). We will find the solution of the Eq. (34) for the boundary conditions of the stable 
states (r — > r or r 2 and r' — > when z — > 00). For these conditions the constant C in the polynom Q(r) is equal to 
0. It is convenient to present the potential Q(r) in the form: 

Q(r) = -Ql(r 2 -r) 2 (r-r a ) 2 , (41) 

where Q\ is the combination of the model parameters. 

The substitution of the expression (41) into the Eq. (34) and integration gives the relation: 



(r-2 - r ) 



In 



r 2 - r(z) 
r - r(z) 



±QeZ. (42) 



After rearrangement of the Eq. (42) we obtain the expression for the internal component of the excitation in the 
following form: 



r(z) = r e + r d th(q e z), 



(43) 



where r e — (r^ + i~q)/2 (r e = r\ for symmetric bistable potential $), = (r% — tq)/2 and q e = ±rdQ e - 

As seen from Eq. (43) at z — > ±00, r{z) — > r2 or ro in dependence on the sing of g e . At z = 0, r(z) = r e (Fig. 4b). 
The width of the excitation is proportional to (r^) -1 . 

For external component in accordance with Eq. (32) and the form of the potential function F(r) in the case of 
equivalent states (Fig. 2b) we have the following: when z — > ±00 then r — > 0, and r(0) = r e = — £-F[r e ] at z = 
(Fig. 4b). Thus the external component has another form in compare with the case of non-equivalent states. In 
the central part of the macromolecule fragment the deformation have the maximum value, and on the edges of the 
excitations the deformation is absent. 

These results may be compared with the data of molecular images of the protein-induced A-B transformation 
of DNA [5]. In the work [5] the structure of the crystal DNA-protein complex [68] was reconstructed in all-atom 
presentations for DNA fragment. As seen from the data of [5], the DNA portion in the complex is in the B-form, 
but free portion of DNA remains in the A-form. At the A-B junction DNA macromolecule bent [5]. The observed 
DNA transformation is in agreement with the results of the present modeling. As in experiment, the static excitation 
for the DNA with equivalent states is caused by the transition between the forms of the double helix, and has the 
maximum of macromolecule deformation at the center of A-B junction (Fig. 4b). 

The obtained results have also a good accordance with the data on intrinsically bent DNA [10,11]. As considered, 
the DNA fragments with the homogenic sequence of AT pair (A-tract) cause the bistability in DNA bending [10- 
12]. In studied DNA the portion with A-tract alter with the DNA portions with the averaged sequence of the base 
pairs. As supposed, the A-tract have two possible conformations: one conformation that is the same as for whole 
macromolecule, and another conformation that differs from the remainder part of double helix. Under some conditions 
the A-tracts transit to another form, and this is considered as the reason of DNA bending [10]. 

The observed DNA bending [10,11] may be explained with the help of the results of the present modeling also. 
According to the case of two equivalent states in the monomer link the bending of DNA fragment is the result of the 
transition from one state to another. 

It is essential that the static excitation in the internal component (43) provides the external component bending 
to one side for both type (there and back) of transitions (Fig. 4b). Therefore, the alternating of the forms in the 
internal component leads to the bending to the side of one groove of the macromolecule. If the length of the excitation 
is divisible to the helix step then the bending of the macromolecule as a whole would realize in one side and in one 
plane. Such mechanism of the intrinsically bending of DNA is in an agreement with the observed properties of the 
bent macromolecule [10-12,69]. 

V. CONCLUSIONS 

The performed study shows the possibility of constructing the relatively simple, two-component models for the 
description of DNA structure transformations on the mesoscopic scales. As was seen for a number of conformational 
transformations in the frame of the double-helical state the model may be formulated in unifying form. It was shown 
that on the pathways of the mesoscopic transformations the DNA macromolecule behaves as homogeneous system in 
the heterogeneous potential field. 

The comparison of the static excitations obtained in the constructed models with the observed deformations of 
the DNA fragments show the qualitative agreement between the theory and experiment. From this correlation 
the conclusion follows that the observed DNA deformations occasioned by intrinsical localized excitations - static 
conformational solitons by their nature. The obtained results may be considered as the direct evidence of the localized 
excitations existence in DNA macromolecule. 
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TABLE I. The parameters of the four-mass model 



Nucleo- 
side 


m* 
(a.u.m.) 


M 4 
(a.u.m.) 


k 

(A) 


e 

(degr.) 


Ade 


203(134) 


312 


5.1 


25 


Thy 


194(125) 


303 


4.6 


35 


Gua 


219(150) 


328 


5.4 


23 


Cyt 


179(110) 


288 


4.5 


30 



* In the bracket there are the masses of the nucleic bases. 



TABLE II. The model parameters for conformational transition in DNA 



Base 


M 


m 






loz 


pair 


(a.u.m.) 


(a.u.m.) 


(a.u.m.) 


(a.u.m.) 


(10- 37 g-cm 2 ) 


AT 


615 


140.7 


99.2 


54.5 


33.7 


GC 


616 


140.8 


98.5 


54.5 


34.5 



TABLE III. The model parameters for DNA opening 



Base 


M 




J 


3 


pair 


(a.u.m.) 


(a.u.m.) 


(10- 37 g-cm 2 ) 


(10" 37 g-cm 2 ) 


AT 


615 


153.7 


5.51 


1.36 


GC 


616 


153.4 


5.80 


1.38 



FIGURE CAPTIONS 

Fig. 1. The four-mass model for the DNA monomer link: (a) - the fragment of the polynucleotide chain and a 
pendulum-nucleoside construction, the nucleoside is shown by the dot-dash line, * - the center of masses of nucleoside; 
(b) - the masses displacements in the plane of the complementary pair; (c) - the double chain of the nucleosides and 
backbone groups masses. 

Fig. 2. A sample view of the potential functions $(r) and F(r): (a) - the case of the double well with non- 
equivalent stable states; (b) - the case of the equivalent states. 

Fig. 3. Static excitation in the case of non-equivalent stable states for DNA monomer link: (a) - the view of the 
effective potential; (b) - the form of the excitations for internal and external components. 



Fig. 4. Static excitation in the case of two equivalent states for DNA monomer link: (a) - the view of the effective 
potential; (b) - the form of the excitations for internal and external components. 



